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Abstract 

We first show that there are in fact triangular arbitrage opportunities in the spot 
foreign exchange markets, analyzing the time dependence of the yen-dollar rate, 
the dollar-euro rate and the yen-euro rate. Next, we propose a model of foreign 
exchange rates with an interaction. The model includes effects of triangular arbitrage 
transactions as an interaction among three rates. The model explains the actual data 
of the multiple foreign exchange rates well. 
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1 Introduction 



The triangular arbitrage is a financial activity that takes advantage of the 
three exchange rates among three currencies. Suppose that we exchange one 
US dollar to some amount of Japanese yen, exchange the amount of Japanese 
yen to some amount of euro, and finally exchange the amount of euro back 
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to US dollar; then how much US dollar do we have? There are opportunities 
that we have more than one US dollar. The triangular arbitrage transaction 
is the trade that takes this type of opportunities. It has been argued that the 
triangular arbitrage makes the product of the three exchange rates converge 
to a certain value [1]. In other words, the triangular arbitrage is a form of 
interaction among currencies. 

The purpose of this paper is to show that there is in fact triangular arbitrage 
opportunities in foreign exchange markets and they generate an interaction 
among foreign exchange rates. We first analyze real data in Sec. 2, showing 
that the product of three foreign exchange rates has a narrow distribution 
with fat tails. In order to explain the behavior, we propose in Sec. 3 a model 
of the time evolution of exchange rates with an interaction. The simulation 
explain the real data well. 



2 Observation of the Triangular Arbitrage Opportunity in High- 
Frequency Data 

In the present paper, we analyze actual data of the yen-dollar rate, the yen- 
euro rate and the dollar-euro rate, taken from January 25 1999 to March 12 
1999 except for weekends. We show in this section that there are actually 
triangular arbitrage opportunities and that the three exchange rates correlate 
strongly. 

2. 1 Existence of Triangular Arbitrage Opportunities 

In order to quantify the triangular arbitrage opportunities, we define the quan- 
tity 

Kt) = f[ n (t), (l) 

i=i 

where ri(t) denotes each exchange rate at time t. We refer to this quantity 
as the rate product. There is a triangular arbitrage opportunity whenever the 
rate product is greater than unity. 

To be more precise, there are two types of the rate product. One is based on 
the arbitrage transaction in the direction of dollar to yen to euro to dollar. 
The other is based on the transaction in the opposite direction of dollar to 
euro to yen to dollar. Since these two values show similar behavior, we focus 
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on the first type of fx(t) hereafter. Thus, we specifically define each exchange 
rate as 



n(t) "yen-dollar ask (t) (2) 

dollar-euro ask (t) ^ 

r 3 (t) = yen-euro bid (t). (4) 

(Note the difference between 

1 [dollar] — > yen-dollar bid [yen] (5) 



and 



1 [yen] yen-dollar ask (6) 



We assume here that an arbitrager can transact instantly at the bid and the 
ask prices provided by information companies and hence we use the prices at 
the same time to calculate the rate product. 

Figure l(a)-(c) shows the actual changes of the three rates: the yen-euro ask, 
the dollar-euro ask and the yen-euro bid. Figure 1(d) shows the behavior of 
the rate product fx(t). We can see that the rate product fx fluctuates around 
the average 

m = (fx(t)) ~ 0.99998. (7) 



(The average is less than unity because of the spread; the spread is the differ- 
ence between the ask and the bid prices and is usually of the order of 0.05% 
of the prices.) The probability density function of the rate product fx (Fig. 2) 
has a sharp peak and fat tails while those of the three rates (Fig. 3) do not. It 
means that the fluctuations of the exchange rates have correlation that makes 
the rate product converge to the average m. 



2.2 Feasibility of the Triangular Arbitrage Transaction 



We discuss here the feasibility of the triangular arbitrage transaction. We 
analyze the duration of the triangular arbitrage opportunities and calculate 
whether an arbitrager can make profit or not. 
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The shaded area in Fig. 2 represents triangular arbitrage opportunities. We 
can see that the rate product is grater than unity for about 6.4% of the time. 
It means that triangular arbitrage opportunities exist about ninety minutes a 
day. The ninety minutes, however, include the cases where the rate product pL 
is greater than unity very briefly. The triangular arbitrage transaction is not 
feasible in these cases. 

In order to quantify the feasibility, we analyze the duration of the triangular 
arbitrage opportunities. Figure 4 shows the cumulative distributions of the 
duration r + of the situation /i > 1 and r_ of /i < 1. It is interesting that the 
distribution of r + shows a power-law behavior while the distribution of r_ dose 
not. This difference may suggest that the triangular arbitrage transaction is 
carried out indeed. 

In order to confirm the feasibility of the triangular arbitrage, we simulate the 
triangular arbitrage transaction using our time series data. We assume that 
it takes an arbitrager T rec [sec] to recognize triangular arbitrage opportunities 
and T ex e[sec] to execute a triangular arbitrage transaction; see Fig. 5. We 
also assume that the arbitrager transacts whenever the arbitrager recognizes 
the opportunities. Figure 6 shows how much profit the arbitrager can make 
from one US dollar (or Japanese yen or euro) in a day. We can see that the 
arbitrager can make profit if it takes the arbitrager a few seconds to recognize 
the triangular arbitrage opportunities and to execute the triangular arbitrage 
transaction. 



3 Modeling 



We here introduce a new model that takes account of the effect of the triangu- 
lar arbitrage transaction as an interaction among the three rates. Many mod- 
els of price change have been introduced so far: for example, the Levy-stable 
non-Gaussian model [2]; the truncated Levy flight [3]; the ARCH/GARCH 
processes [4,5]. They discuss, however, only the change of one price. They did 
not consider an interaction among multiple prices. As we discussed in Sec. 
2, however, the triangular arbitrage opportunity exists in the market and is 
presumed to affect price fluctuations in the way the rate product tends to 
converge to a certain value. 
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3.1 Basic Time Evolution 



The basic equation of our model is a time-evolution equation of the logarithm 
of each rate: 

\n n {t + At) = lnr,(t) + f t (t) + g{y(t)), (i = 1, 2, 3) (8) 

where v is the logarithm of the rate product 

u(t) = lnf,(t) = j2lnr t (t). (9) 

i=i 

Just as /x fluctuates around m — (/j) ~ 0.99998, the logarithm rate product z/ 
fluctuates around 

e = (In y) ~ -0.00091 (10) 



(Fig. 7(a)). In this model, we focus on the logarithm of the rate-change ratio 
ln(rj(t + At)/ri(t)), because the relative change is presumably more essen- 
tial than the absolute change. We assumed in Eq. (8) that the change of the 
logarithm of each rate is given by an independent fluctuation fi(t) and an 
attractive interaction g{v). The triangular arbitrage is presumed to make the 
logarithm rate product v converge to the average e; thus, the interaction func- 
tion g(y) should be negative for v greater than e and positive for v less than 
e: 



, < , for v > e 

9(") { (11) 
> , for v < e. 

As a linear approximation, we define g(y) as 

g(y) = - a (y- e ) ( i2 ) 

where a is a positive constant which specifies the interaction strength. 
The time-evolution equation of v is given by summing Eq. (8) over all v. 

v{t + At) - e = (1 - 3a)(is(t) - e) + F(t), (13) 
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where 



^(0 = E /<(*)■ (14) 

1=1 



This is our basic time-evolution equation of the logarithm rate product. 
3.2 Estimation of Parameters 



The interaction strength a is related to the auto-correlation function of v as 
follows: 

1 - 3a = c(At) = (z/2(t)) _ (z/(t))2 • (15) 



Using Eq. (15), we can estimate a(At) from the real data series as a function 
of the time step At. The auto-correlation function c(At) is shown in Fig. 8(a). 
The estimate of a(At) is shown in Fig. 8(b). Hereafter, we fix the time step 
at At = 60 [sec] and hence use 

a(l[min]) = 0.17 ±0.02 (16) 



for our simulation. 

On the other hand, the fluctuation of foreign exchange rates is known to be a 
fat-tail noise [6,7]. Here we take fi(t) as the truncated Levy process [3,8]: 

Ptlf(/;«, 7 ,0 = qP L (f;a n )Q(l- |/|), (17) 

(18) 

where q is the normalization constant, Q(x) represents the step function and 
Pl(x; a, 7) is the symmetric Levy distribution of index a and scale factor 7: 

1 °° 

P h {x- a, 7) = - / e-^ a cos(kx)dk < a < 2. (19) 

7T J 





We determine the parameters a, 7 and / by using the following relations for 
1 < a < 2 [6,9]: 



a(a — 1)7 , . , 

ca= | S ; 9 ^ 2 - a , 20 
I cos(7TQ;/2)| 
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rate 



a 



7 



n (1 /yen-dollar ask) 1.8 7.61xlCT 7 1.38xlCT 2 
r 2 (1/dollar-euro ask) 1.7 4.06xl(T 7 3.81xl(T 2 
r 3 (yen-euro bid) 1.8 6.97xl(T 7 7.58xl(T 2 



Table 1 

The estimates of the parameters. 

_ (3-a)(2-a)|cos(7ra/2)| 



(21) 




where c n denotes the nth cumulant and k is the kurtosis k = c 4 /c 2 2 . The 
estimates are shown in Table 1. The generated noises with the estimated 
parameters are compared to the actual data in Fig. 9. 

We simulated the time evolution (13) with the parameters given in Eqs. (10), 
(16) and Table 1. The probability density function of the results (Fig. 7(b)) 
is compared to that of the real data (Fig. 7(a)) with At = l[min] in Fig. 
10. The fluctuation of the simulation data is consistent with that of the real 
data. In particular, we see good agreement around v ~ e as a result of the 
linear approximation of the interaction function. Figure 7(c) shows u{t) of the 
simulation without the interaction, i.e. a = 0. The quantity v fluctuates freely, 
which is inconsistent with the real data. 



3.3 Analytical approach 

We can solve the time-evolution equation (13) analytically in some cases. Let 
us define 



u{t) - e 



(22) 



and 



< 



F 



c 



u>(t), u>' 
F(t), F' 
c(At). 



u(t + At), 
F(t + At), 



(23) 



Equation (13) is then reduced to 



u' = cu + F 



(24) 
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Assume that the probability of uj having a value in uj ~ uj + duj is P u ((v). The 
joint probability of uj' having a value in uj' ~ uj' + duj' and F' having a value 
in F' ~ F' + dF' is given by 

P^,f'(oj', F')duj'dF' = P^PpiF^dujdF' 

= -PU^—^)PF(F')<kj'dF', (25) 
c c 

where Pf{F) is the probability of F having a value in F ~ F + dF. The 
probability density function of uj' is thus given by 

PA«S) = ~J PU^^)P F (F')dF'. (26) 

The characteristic function of uj' is the Fourier transform 

PAk) = J PA^)^' k ^' = Pu,(ck)P F (k), (27) 

where P w and P F are the Fourier transforms of P w and P F) respectively. Then 
we obtain 

Pw(t) (k) = P u (t-At) (ck)P F ( t -At) (k) 

N 

= P^ 0) (c N k)l[P F (c n k) 

n=0 
3 N 

= P^)(c N k)UUPfAc n k) 

i=l n=0 

3 N 

= UUh(c n k), (28) 

1=1 71=0 

where t = AT At. We here assumed P^q^lj) = 5(uj) and hence P u (o)(k) = 1. 

The above argument shows the essential reason of the sharp peak and fat tails 
in Fig. 2. If we had c = 1, or a = (without the interaction), the noise F at 
every time step would accumulate in uj and the probability density function 
of uj = v — e would be Gaussian due to the central limit theorem. If we have 
c < 1, or a > (with the interaction), the noise at the past time steps decay 
as c n . The largest contribution to uj comes from the noise one time step before, 
which is a fat-tail noise [3] . 

As a special case, if the noises fi(t) obey a Levy distribution of the same index 
a and the same scale factor 7, namely if 

P h (k) = e" 7|fc|Q foralH, (29) 



8 



the distribution of & is also a Levy distribution of the same index a and a 
different scale factor 7' given by 



7 ' l-{c(At)}« T 



4 Summary 

We first showed that triangular arbitrage opportunities exist in the foreign 
exchange market. The rate product \i fluctuates around a value m. Next, 
we introduced a model including the interaction caused by the triangular 
arbitrage transaction. We showed that the interaction is the reason of the 
sharp peak and the fat-tail property of the distribution of the logarithm rate 
product v. Finally we showed that our model is solvable analytically in some 
cases. 
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Fig. 1. The time dependence of (a) the yen-dollar ask l/r±, (b) the dollar-euro ask 
l/ r 2, (c) the yen-euro bid r% and (d) the rate product [i. The horizontal axis denotes 
the seconds from 00:00:00, January 25 1999. 
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(a) 



(b) 
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Fig. 2. The probability density function of the rate product ji. (b) is a 
semi- logarithmic plot of (a). The shaded area represents triangular arbitrage op- 
portunities. The data were taken from January 25 1999 to March 12 1999. 
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Fig. 3. The probability density function of the three rates: (a) the reciprocal of the 
yen-dollar ask, n, (b) the reciprocal of the dollar-euro ask, r2 and (c) the yen-euro 
bid T3. The data were taken from January 25 1999 to March 12 1999. 
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Fig. 4. The cumulative distributions of r + (o) and r_ (o). The distribution of r + 
shows a power-low behavior. The data were taken from January 25 1999 to March 
12 1999. 
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Fig. 5. A conceptual figure of the profit calculation. We assume that it takes an 
arbitrager T rec [sec] to recognize triangular arbitrage opportunities and T exe [sec] to 
execute a triangular arbitrage transaction. The circles (o) indicate the instances 
where triangular arbitrage transactions are carried out. 




Fig. 6. A phase diagram of the profit that an arbitrager can make from one US 
dollar (or Japanese yen or euro) in a day under the assumption shown in Fig. 5. 
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Fig. 7. The time dependence of ^(t[min]) of (a) the real data, (b) the simulation data 
with the interaction and (c) without the interaction. In (b), v fluctuates around e 
like the real data. 
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Fig. 8. (a) The auto-correlation function of u, c(At). (b) The time-step dependence 
of o(At). 
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Fig. 9. The cumulative distributions of the one-minute changes 
|lnrj(t + l[min]) — lnrj(t)| (o represents upward movements and o repre- 
sents downward movements) and the generated noise /j ( — ): (a) the yen-dollar ask 
and fi, (b) the dollar-euro ask and /2, and (c) yen-euro bid and ^3. The real data 
were taken from January 25 1999 to March 12 1999. 
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Fig. 10. The probability density function of v. The circle (o) denotes the real data 
and the solid line denotes our simulation data with the interaction. The simulation 
data fit the real data well. 
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